Reduced-order modeling for unsteady transonic flows around an airfoil

High-transonic unsteady flows around an airfoil at zero angle of incidence and moderate Reynolds numbers are characterized by an unsteadiness induced by the von Kármán instability and buffet phenomenon interaction. These flows are investigated by means of low-dimensional modeling approaches. Reduced-order dynamical systems based on proper orthogonal decomposition are derived from a Galerkin projection of two-dimensional compressible Navier-Stokes equations. A specific formulation concerning density and pressure is considered. Reduced-order modeling accurately predicts unsteady transonic phenomena

Reduced-order modeling of transonic flows around an airfoil submitted to small deformations

A reduced-order model (ROM) is developed for the prediction of unsteady transonic flows past an airfoil submitted to small deformations, at moderate Reynolds number. Considering a suitable state formulation as well as a consistent inner product, the Galerkin projection of the compressible flow Navier–Stokes equations, the high-fidelity (HF) model, onto a low-dimensional basis determined by Proper Orthogonal Decomposition (POD), leads to a polynomial quadratic ODE system relevant to the prediction of main flow features. A fictitious domain deformation technique is yielded by the Hadamard formulation of HF model and validated at HF level. This approach captures airfoil profile deformation by a modification of the boundary conditions whereas the spatial domain remains unchanged. A mixed POD gathering information from snapshot series associated with several airfoil profiles can be defined. The temporal coefficients in POD expansion are shape-dependent while spatial POD modes are not. In the ROM, airfoil deformation is introduced by a steady forcing term. ROM reliability towards airfoil deformation is demonstrated for the prediction of HF-resolved as well as unknown intermediate configurations

Anisotropic Organised Eddy Simulation for the prediction of non-equilibrium turbulent flows around bodies

The unsteady turbulent flow around bodies at high Reynolds number is predicted by an anisotropic eddy-viscosity model in the context of the Organised Eddy Simulation (OES). A tensorial eddy-viscosity concept is developed to reinforce turbulent stress anisotropy, that is a crucial characteristic of non-equilibrium turbulence in the near-region. The theoretical aspects of the modelling are investigated by means of a phase-averaged PIV in the flow around a circular cylinder at Reynolds number 1.4×10^5. A pronounced stress–strain misalignment is quantified in the near-wake region of the detached flow, that is well captured by a tensorial eddy-viscosity concept. This is achieved by modelling the turbulence stress anisotropy tensor by its projection onto the principal matrices of the strain-rate tensor. Additional transport equations for the projection coefficients are derived from a second-order moment closure scheme. The modification of the turbulence length scale yielded by OES is used in the Detached Eddy Simulation hybrid approach. The detached turbulent flows around a NACA0012 airfoil (2-D) and a circular cylinder (3-D) are studied at Reynolds numbers 105 and 1.4×10^5, respectively. The results compared to experimental ones emphasise the predictive capabilities of the OES approach concerning the flow physics capture for turbulent unsteady flows around bodies at high Reynolds numbers

Deflection, drift and advective growth in variable-density, laminar mixing layers

Specific features of the variable-density mixing layers without gravity effects are studied using self-similar solutions to the laminar and time-evolving variant of this flow. Density variations come from either mass or temperature mixing, accounting for, in the latter case, the effect of the Mach number. The transverse profiles of the flow quantities, as well as the time evolutions of the global characteristic scales of the mixing layer, are given for a wide range of density ratio and Mach-number values. When compared to the constant-density case, it appears that most of the specificity of these flows comes from the emergence of a nonzero transverse component of the velocity. First, it produces a deflection of the flow that can be either confined in the core of the layer or global, the whole layer being tilted at an angle from the initial flow direction. In most cases, this deflection is such that some part of the higher-density fluid is "entrained" in the direction of the lower-density fluid, leaving no possibility to define a dividing streamline. Second, it leads to a shift between the density profile and the profiles of the other flow quantities. This shift scales on the time-increasing mixing-layer thickness and therefore appears as a time drift. When global deflection is present, the tilting of the layer can be shown to be equivalent to a global drift of the mixing/shear layer toward the light-fluid side of the flow. Third, transport by the transverse velocity component affects the spreading of the mixing layer, giving rise to an additional effect referred to as advective growth. Examination of the density-ratio and Mach-number effects leads to surprising results: While the momentum thickness is always observed to decrease when increasing these parameters, conventional thicknesses based on the profiles of the different variables can show opposite behaviors depending on the form of the diffusion model for the considered variable

Three-dimensional flow past a fixed or freely vibrating cylinder in the early turbulent regime

The three-dimensional structure of the flow downstream of a circular cylinder, either fixed or subjected to vortex-induced vibrations, is investigated by means of numerical simulation, at Reynolds number 3900, based on the cylinder diameter and current velocity. The flow exhibits pronounced fluctuations distributed along the span in all studied cases. Qualitatively, it is characterized by spanwise undulations of the shear layers separating from the body and the development of vortices elongated in the plane normal to its axis (planar vortices). A quantitative analysis of crossflow vorticity fluctuations in the spanwise direction reveals a peak of fluctuation amplitude in the near region (i.e., area of formation of the spanwise wake vortices) and opposite trends of the spanwise wavelength in the shear layer and wake regions; the wavelength tends to decrease as a function of the streamwise distance in the shear layers down to a minimum value close to 0.5 body diameters and then slowly increases further in the wake. The spanwise structure of the flow is differently altered in these two regions, once the cylinder vibrates. In the shear layer region, body motion is associated with an enhancement of planar vortex formation. The amplification of vorticity spanwise fluctuations in this region is accompanied by a reduction of the spanwise wavelength; it is found to decrease as a function of the instantaneous Reynolds number based on the instantaneous flow velocity seen by the moving body, following the global trend of the wavelength versus Reynolds number previously reported for fixed cylinders. In the wake region, the flow spanwise structure is essentially unaltered compared to the fixed body case, in spite of the major distortions of the streamwise and crossflow length scales

Vortex-induced vibrations of a cylinder in planar shear flow

The system composed of a circular cylinder, either fixed or elastically mounted, and immersed in a current linearly sheared in the cross-flow direction, is investigated via numerical simulations. The impact of the shear and associated symmetry breaking is explored over wide ranges of values of the shear parameter (non-dimensional inflow velocity gradient, β ∈ [0, 0.4]) and reduced velocity (inverse of the non-dimensional natural frequency of the oscillator,U∗ ∈ [2, 14]), at Reynolds number Re = 100; β, U∗ and Re are based on the inflow velocity at the centre of the body and on its diameter. In the absence of large-amplitude vibrations and in the fixed body case, three successive regimes are identified. Two unsteady flow regimes develop for β ∈ [0, 0.2] (regime L) and β ∈ [0.2, 0.3] (regime H). They differ by the relative influence of the shear, which is found to be limited in regime L. In contrast, the shear leads to a major reconfiguration of the wake (e.g. asymmetric pattern, lower vortex shedding frequency, synchronized oscillation of the saddle point) and a substantial alteration of the fluid forcing in regime H. A steady flow regime (S), characterized by a triangular wake pattern, is uncovered for r β > 0.3. Free vibrations of large amplitudes arise in a region of the parameter space that encompasses the entire range of β and a range of U∗ that widens as β increases; therefore vibrations appear beyond the limit of steady flow in the fixed body case (β = 0.3). Three distinct regimes of the flow-structure system are encountered in this region. In all regimes, body motion and flow unsteadiness are synchronized (lock-in condition). For β ∈ [0, 0.2], in regime VL, the system behaviour remains close to that observed in uniform current. The main impact of the shear concerns the amplification of the in-line response and the transition from figure-eight to ellipsoidal orbits. For β ∈ [0.2, 0.4], the system exhibits two well-defined regimes: VH1 and VH2 in the lower and higher ranges of U*, respectively. Even if the wake patterns, close to the asymmetric pattern observed in regime H, are comparable in both regimes, the properties of the vibrations and fluid forces clearly depart. The responses differ by their spectral contents, i.e. sinusoidal versus multi-harmonic, and their amplitudes are much larger in regime VH1, where the in-line responses reach 2 diameters (0.03 diameters in uniform flow) and the cross-flow responses 1.3 diameters. Aperiodic, intermittent oscillations are found to occur in the transition region between regimes VH1 and VH2; it appears that wake-body synchronization persists in this case

Capturing transition features around a wing by reduced-order modeling based on compressible Navier-Stokes equations

The three-dimensional transition in the flow around a NACA0012 wing of constant spanwise section at Mach number 0.3, Reynolds number 800, and incidence 20° is investigated by direct numerical simulation and reduced-order modeling. The interaction between the von Kármán and the secondary instabilities is analyzed. Irregular events in the flow transition modulating the spanwise undulation are highlighted and quantified. These transition features, including "local intermittencies" in the secondary instability pattern, are efficiently captured by a reduced-order model derived by means of the Galerkin projection of the compressible flow Navier-Stokes equations onto a truncated proper orthogonal decomposition basis

On the multidisciplinary control and sensing of a smart hybrid morphing wing

Morphing wing technology is of great interest for improving the aerodynamic performance of future aircraft. A morphing wing prototype using both surface embedded Shape Memory Alloys (SMA) and piezoelectric macro fiber composite (MFC) actuators has been designed for wind tunnel experiments. This smart wing is a mechatronic system that contains embedded sensors to measure the surrounding flow and control the actuators. This article will focus on the control of the cambering system which is achieved using a group of nested control loops as well as on the perspective of a novel control strategy using in-situ temperature measurements. It will be shown that by exploiting the inherent hysteretic properties of the SMAs cambering a significant reduction in power consumption is possible by appropriately tailoring the control strategy. Furthermore, by comparing the post-processed pressure signals recorded during the wind tunnel experiments to the aerodynamic performance gains a perspective for a novel in-situ control will be shown

Simulation numérique et analyse physique du tremblement transsonique d'un profil supercritique à Reynolds élevé

L'écoulement autour d'un profil d'aile supercritique est étudié numériquement dans le régime du tremblement transsonique à Reynolds élevé par des approches de modélisation de la turbulence statistiques (URANS), statistiques avancées (Organized-Eddy Simulation) et hybrides (DDES). Les différentes méthodes sont évaluées par rapport à leurs capacités de prédiction des caractéristiques du mouvement de l'onde de choc, des fluctuations de grandeurs locales de l'écoulement et des profils de vitesse sur l'extrados. L'influence de l'instabilité de von Kármán sur le tremblement est également analysée

Two-degree-of-freedom vortex-induced vibrations of a circular cylinder at Re=3900

The vortex-induced vibrations of an elastically mounted circular cylinder are investigated on the basis of direct numerical simulations. The body is free to move in the in-line and cross-flow directions. The natural frequencies of the oscillator are the same in both directions. The Reynolds number, based on the free stream velocity and cylinder diameter, is set to 3900 and kept constant in all simulations. The behavior of the coupled flow-structure system is analyzed over a wide range of the reduced velocity (inverse of the natural frequency) encompassing the lock-in range, i.e. where body motion and flow unsteadiness are synchronized. The statistics of the structural responses and forces are in agreement with prior experimental results. Large-amplitude vibrations develop in both directions. The in-line and cross-flow oscillations are close to harmonic; they exhibit a frequency ratio of 2 and a variable phase difference across the lock-in range. Distinct trends are noted in the force-displacement phasing mechanisms in the two directions: a phase difference jump associated with a sign change of the effective added mass and a vibration frequency crossing the natural frequency is observed in the cross-flow direction, while no phase difference jump occurs in the in-line direction. Higher harmonic components arise in the force spectra; their contributions become predominant when the cylinder oscillates close to the natural frequency. The force higher harmonics are found to impact the transfer of energy between the flow and the moving body, in particular, by causing the emergence of new harmonics in the energy transfer spectrum